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Propagation  of  Hydromagnetlc  Waves 

The  general  problem  Is  the  description  of  the  properties 
of  small  amplitude  signals  propagating  in  a  plesma.   To  this 
end,  the  computation  of  a  dispersion  relation  should  bo  con- 
sidered to  be  the  beginning  rather  than  the  end  of  the 
problem.   In  specific  cases  which  could  be  evaluated  in 
detail,  the  properties  of  plane  waves  were  found  to  be  not 
only  incomplete  but  even  misleading.   For  example,  non-plane 
fronts  can  be  found  which  propagate  faster  in  some  directions 
than  the  maximum  speed  of  propagation  of  a  plane  wave  in  any 
direction  [22-VIII];  similarly,  one  mode  whose  plane  wave 
fronts  propagate  slower  than  those  of  another  mode  in  every 
direction  can  catch  up  and  pass  the  second  mode  if  the  front 
is  not  plane  [33.   Complete  modes  (as  distinguished  from 
plane  waves)  are  found  to  separate  only  in  special  cases  and, 
even  then,  only  in  a  certain  subtle  sense  [31. 

The  types  of  question  which  were  approached  are  the 
possibility  of  separation  of  modes  (in  arbitrary  disturbances 
rather  than  plane  waves),  channeling  of  energy  propagation 
and  radiation  patterns,  and  the  general  propagation  of  arbitrary 
wave  fronts  including  refraction,  focussing,  etc,  with  the 
idea  of  eventual  application  to  localized  plasma  diagnostics. 
The  deeper  questions  force  one  to  study,  at  least  initially, 
some  of  the  simpler  models  of  plasma  behavior,   Eor  the 
classical  perfectly  conducting  fluid,  an  extensive  and  almost 
complete  analysis  has  been  made  of  the  propagation  of  arbitrary 
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wave  fronts  Into  a  general  inhomogeneous  medluin  (geometrical 
optics) [18 j.   Not  only  the  position  of  the  front,  but  also 
the  variation  in  amplitude  are  given  by  simple  formulas.   The 
more  difficult  problem  of  general  energy  and  signal  transfer 
(not  necessarily  related  to  a  wave  front)  has  been  success- 
fully attacked  in  special  geometries  [3]  and  is  being  extended 
to  more  general  cases  by  the  use  of  recently  refined  asymptotic 
techniques  119].   For  perturbations  about  a  constant  magnetic 
field,  one  mode  was  found  which  propagates  purely  one- 
diraensionally,  guided  along  the  magnetic  lines  [3].   Con- 
ceivably, this  could  be  used  as  a  diagnostic  tool  to  obtain 
average  properties  of  an  open-ended  plasma  along  a  single 
magnetic  line.   The  oroblem  of  propagation  into  a  non-iiniform 
field  is  being  studied. 

Extension  of  these  results  to  more  realistic  plasma 
models  using  the  Boltzmann  equation  are  \inder  way.   Since 
many  of  the  questions  have  never  been  satisfactorily  answered 
for  even  a  simple  gas  without  ionization  or  electromagnetic 
effects  present,  it  was  found  to  be  necessary  to  first  develop 
some  parts  of  this  simpler  theory  before  dealing  with  plasmas. 
As  a  first  step,  an  asymptotic  approximation  was  obtained  to 
the  time-independent  fundamental  solution  of  the  Boltzmann 
equation  [Ij.],   This  is  being  extended  to  time-dependent 
propagation  and  to  similar  analysis  in  plasmas. 
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Particle  Orbits 

The  proof  of  the  asymptotic  validity  of  the  guiding 
center  expansion  [2?]  was  submitted  for  publication  [10], 

A  large  number  of  general  results  concerning  adiabatic 
invariants  and  the  asymptotic  behavior  of  charged  particle 
orbits  was  obtained.   In  September  (this  was  annoiinced 
informally  at  the  Geneva  conference)  C,  S.  Gardner  proved 
the  existence  of  the  "second"  adiabatic  invariant  (even  to 
arbitrarily  high  order)  for  axially  symmetric  systems.   These 
results  were  extended  [Ij-Oa],  using  classical  perturbation 
theory  of  Hamlltonian  systems  to  obtain 

1)  a  new  and  simplified  proof  of  the  invarlance  of  the 
magnetic  moment  to  arbitrarily  high  order 

2)  under  appropriate  circumstances  (roughly  some  sort  of 
mirror  arrangement  --  but  not  necessarily  a_xlally  symmetric) 
the  invariance  to  arbitrary  order  of  a  second  adiabatic 
invariant. 

3)  \''aider  appropriate  circumstances  the  invarlance  of  a 
third  invariant 

I4.)   extension  of  all  these  results  to  relativistic 
particles. 

The  Hamlltonian  treatment  yields  a  reduction  of  the  original 
sixth  order  system  to  fourth  order  on  elimination  of  the 
magnetic  moment  as  a  "constant"  of  the  motion  and  a  similar 
reduction  to  second  order  on  elimination  of  the  second  invar- 
iant.  What  remains  is  a  system  of  one  degree  of  freedom 
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describing  the  motion  of  the  "line"  to  which  the  particle  is 
"tied".  A  brief  accoiint  was  submitted  for  publication  [13], 
and  a  more  detailed  version  is  being  prepared  for  publication. 

These  adiabatic  results  do  not  give  much  information 
about  deviations  from  adiabaticity.   Some  preliminary  analytic 
results  in  this  direction  were  obtained,  but  much  more  work 
is  required. 

Special  problems  in  which  aopllcations  of  particle  orbit 
theory  arise  are  described  under  Shock  Theory,  Equilibrium  and 
Stability,  and  Cusped  Geometries. 


Shock  Theory 

The  analysis  of  the  self-consistent  steady  shock  (non- 
colliding  particles)  propagating  in  a  direction  perpendicular 
to  the  magnetic  field  was  continued.   It  had  previously  been 
observed  that  there  are  at  least  two  significantly  different 
sources  of  entropy  production,  in  particular  one  originating 
from  particles  which  are  reflected  back  by  their  first 
encounter  with  the  charge  separation  electric  field  [ic].   It 
was  now  observed  that  there  are  parameter  values  for  which 
this  entropy  increase  is  an  appreciable  part  of  the  maximum 
which  could  occur  across  a  complete  shock  [1+Oc].   Thus,  even 
ignoring  other  mechanisms,  there  can  be  appreciable  irrevers- 
ible heating  of  the  ions  in  a  very  thin  layer.   Specifically, 
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there  are  interesting  ranges  of  the  physical  parameters 
which  yield  ar  essentially  discontinuous  shock  front  across 
which  ions  alone  are  irreversibly  heated;  this  front  is 
followed  by  an  approximately  adiabatic  periodic  wave  train. 
A  full-scale  analysis  and  computation  of  the  magnitude  of 
this  heating  as  a  function  of  various  parameters  vjas  carried 
out  and  will  be  Issued  shortly  [11]. 

Details  of  some  of  the  earlier  work  in  this  subject  were 
Issued  in  [2], 

Several  fluid-type  approximations  to  this  problem  were 
continued  (cf,  [35g]),  but,  thus  far,  unsuccessfully  (i.e.,  in 
prediction  of  a  damped  shock  profile). 


Transport  Theory 

By  making  strong  use  of  symmetry  relations,  it  was 
found  to  be  possible  to  elucidate  and  clarify  the  jungle  of 
transport  coefficients  and  collision  integrals  which  arise 
from  the  Boltzmann  equation  [1L|.].   The  program  aimed  at  im- 
proved values  for  these  coefficients,  see  [29]  and  [?],  was 
continued.   A  dictionary  of  collision  integrals  is  being 
prepared  (L.  Sirovich)  to  avoid  future  repetitions  of  the 
extensive  redundant  calculations  in  the  literature. 
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The  question  of  the  slovjing  down  and  thermalization  of 
a  single  fast  ion  was  treated,  [l^Odj,  [12].   Care  had  to  be 
taken  to  include  both  electron  and  ion  effects  since  either 
one  dominated  in  different  regions.   Quite  complicated 
thermalization  patterns  were  found.  Including  development  of 
a  secondary  peak  in  the  distribution  and  the  transition  from 
an  initial  6-function  through  a  distribution  much  wider  than 
the  ultimate  Maxwellian  before  final  thermalization.   The 
final  decav  to  Maxwellian  was  found  to  be  approximately 
exponential,  but  the  decay  time  depends  strongly  on  the  initial 
energy  of  the  ion.   In  other  words,  reference  to  a  universal 
ion  "thermalization"  time  can  be  quite  misleading. 

The  asymptotic  evaluation  of  the  tail  end  of  the  ion 
distribution  and  resultant  correction  In  the  thermonuclear 
reaction  rate  [28],  [le],  was  improved  to  obtain  an  exact 
rather  than  an  asymptotic  result  (R.  Liboff  -  nixmerical  results 
not  yet  complete). 

Two  questions  of  anomalous  resistivity  induced  by  plasma 
oscillations  were  considered.   The  first  depends  on  the  increased 
dispersion  of  a  particle  due  to  the  presence  of  externally 
induced  plasma  oscillations  which  are  only  slowly  damped  [L|.Ob]. 
This  effect  can  be  exceptionally  large  in  a  sheath.   The  second 
results  from  the  introduction  of  plasma  oscillations  into  the 
formulation  of  Ohm's  law,  thereby  modifying  the  diffusion  rate 
of  plasma  into  field  (J.  Thurber,  not  yet  complete).   These 


-  9  - 


two  effects  are  independent;  in  particular,  the  plasma  resis- 
tivity value  to  be  used  in  the  second  problem  should  be  the 
corrected  value  as  given  by  the  first. 


Equilibrium  and  Stability 

The  problem  of  the  self-consistent  equilibrium  of  a 
mirror-machine  geometry  in  the  guiding  center  limit  (with 
non-diagonal  stress  tensor,  generalizing  the  fluid  model, 
Vp  =  J  X  B)  was  solved  in  certain  limiting  cases  (A,  Oppenheim 
report  to  appear).   More  general  solutions  also  seem  to  be 
accessible  and  are  being  studied.  A  stability  analysis  will 
follow. 

Existence  theorems  have  been  proved  for  a  variety  of 
axially  symmetric  and  toroidal  free  boundary  equilibria  [l^Oe] 
(M.Schechter  and  C.  Morawetz  -  detailed  report  to  appear). 

The  problem  of  the  stability  of  a  toroidal  "stabilized" 
pinch  has  been  formulated  and  prepared  for  computation  (R,  D, 
Richtmyer,  R.  Lust,  B.  Suydam). 

A  stability  analysis  using  the  guiding  center  fluid 
equations  (with  non-diagonal  stress  tensor)  [3i4ii]  was  found 
to  yield  instability  when  the  ratio  of  parallel  to  perpen- 
dicular pressure  is  either  too  large  or  too  small  [l4.0f].   At 
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second  glance,  this  result  is  found  to  be  illusory  or,  at 
least,  hard  to  interpret  since  the  ecuations  of  motion  are 
found  to  be  not  well-posed  when  the  stated  inequalities  are 
exceeded.   In  other  words,  the  motion  is  not  unstable;  it  is 
nonexistent.   This  result  may  cast  doubt  on  all  stability 
analyses  which  use  the  guiding  center  limit  since  it  is  this 
limit  which  seems  to  produce  the  paradoxical  result. 

Comprehensive  and  detailed  reports  covering  the  material 
which  was  briefly  summarized  in  the  two  Geneva  papers  [lb] 
and  [Id]  will  be  issued  shortly,  [8-IX]  and  [8-X]. 


Cusped  Geometries 

Work  In  this  subject  consisted  mainly  in  the  continuation 
of  problems  which  have  been  under  study  for  some  time.   These 
include  the  computation  of  axially  symmetric  free  boxondary 
configurations  and  the  computation  of  nonadiabatic  orbits  in 
cusped  magnetic  fields  for  the  improved  estimation  of  particle 
loss  rates. 

References:   [la],  [25],  [6],  [9]. 


Experimental 

The  experimental  program  was  terminated  as  of  July  19^9. 
A  final  report  is  being  prepared  [151. 
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